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We consider the stochastic ranking process with the jump times of the particles 
determined by Poisson random measures. We prove that the joint empirical dis- 



•n , tribution of scaled position and intensity measure converges almost surely in the 

. infinite particle limit. We give an explicit formula for the limit distribution and 



show that the limit distribution function is a unique global classical solution to 
an initial value problem for a system of a first order non-linear partial differential 
O . equations with time dependent coefficients. 



fXj ■ 1 Introduction. 



Let A^(M+) be the space of Radon measures p on the Borel cx-algebra i3(]R+) of non- 
negative reals M+. Let be a positive integer, and let z^f^'', i = 1, 2, . . . , A^, be indepen- 
dent Poisson random measures (Poisson point processes) on ]R_|_ , defined on a probability 
space (P, J-", fi). For each i, denote the intensity measure of u^^^ by pf^^; 

(1) E[4''\A)]=pf\A), AeB{R^). 
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Throughout this paper, we assume pf^'^ G A^(]R+) and that pf^ is continuous (i.e., 
p(^)({t}) = 0, t ^ 0) for all N and i. 

Let x{^\x^2^\ . . . , x'^^^ be a permutation of 1, 2, . . . , A^, and define a process 



xw = (xr,...,xf)^ 



"AT 



by 
(2) 



k=i Jo ' Jo 



i = l,2,...,N, t>0. 



where, 1a is the indicator function of event A. 

Denote the unit measure concentrated on c by 5c . With probability 1 we can write 

oo 

(3) z.f)=^5(,), ^ = 1,2,...,N, 

where, with probability 1, Tiff's are random variables satisfying < t-^^ < t-^^ < ■ ■ • , 
i = 1,2, . . . , N , and Tj;^^ ^ t-,^J if 7^ [i' ■,]')■ In the following, we work on the event 
that these inequalities hold. 

The right hand side of ([2]) is a simple function in t. At t = Xj ,,■ we see 

.ff'(r„)-A-."(r„-) = l-Xf'(r.,-), 

which implies 

(4) xf\n,) = 1. 

With similar consideration, we see that the process X^^^ is uniquely determined by 
Explicitly, we have, for i = 1, . . . , N, 

SN) 



(5) xr\t) 



TV 

^ + Yl ^'eN; r.(^'<rW<i = ^ < J = 1, 2, 3, ... . 



In the case of the (homogeneous) Poisson process (i.e., the case p-^-*((0,t]) = wf^h, 
t ^ 0, for positive constants w\^^), a discrete time version of the process ([5]) has been 
known for a long time [251 ISSl [IH [6l [211 [IH] and is called move-to-front (MTF) rules. 
The process has, in particular, been extensively studied as a model of least-recently-used 
(LRU) caching in the field of information theory [231[Sl[a[Zl[Sl[2a[SllIIllini[ni 
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and also is noted as a time-reversed process of top-to-random shuffling. With a great 
advance in the internet technologies, a new application of the process appeared [15] . 
The ranking numbers such as those found in the web pages of online bookstores are found 
to follow the predictions of the model. 

In [12], the case where ^'^^'''s are (homogeneous) Poisson processes with pl^\{0,t]) = 
w^^^t is considered, and the joint empirical distribution of jump rate w^'^^ and normalized 
position 

(6) Yrit) = ^iXr\t)-l), 
1 ^ 

given by /ij^"* = y(^)(t)) ? is studied. (We will abuse notation slightly and 

i=l 

denote a unit measure on any space by 6c.) It is proved in [12] that a scaling limit 

(7) fit = lim fif^ 

N^oo 

exists (under reasonable assumptions), and an explicit formula for /if, which is a deter- 
ministic distribution on R_|_ x [0, 1), is given. In [13], it is proved that, if the scaling limit 
of the jump rate distribution is a discrete distribution, the limit fit is the unique time 
global solution to an initial value problem for a system of first order non-linear partial 
differential equations (inviscid Burgers equations with a term representing evaporation). 
The structure of the explicit formula for fit is naturally explained by a standard method 
of characteristic curves for the solution to the partial differential equations. 

In the present paper, we will generalize the main results of [121 113] to the case where 
z/|^'''s are Poisson random measures. We shall call the process X'^'^^ defined by ([2]), or 
equivalently by ([5]), a stochastic ranking processes after [T2l [T3l [T4] . 

Put 

N 

(8) X^c\t) = T.K^-^<v t^^- 

' i.l — 

i=l 

XQ^\t) is a random variable which denotes the position of the boundary between the 
top side X ^ X^\t) and the tail side x > X^^t), where each particle in the top side 
(i.e., i which satisfies X^^\t) ^ xjj^\t)) has experienced jump to the top by time t (i.e., 
t^'t^^ ^ t), and the particles in the tail side are those particles which have not jumped to 
the top by time t. 

Proposition 1.1. Lett ^ 0, and assume that a sequence of distributions {Xi^^ ; N G 
N} on ]R_|_ defined by 

1 ^ 

1=1 
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converges weakly as N oo to a probability distribution At . Then the scaled position of 
the boundary 



N 



(10) i^"'w = ^-4"'w = ^Ei.<-s. 

converges almost surely as N oo to 

(11) yc{t) = 1-1 e-'Xt{ds). 







O 

Proof. The definition implies that Ylf\t) - E[ Y^f^\t) ] is an arithmetic mean 
of independent variables 



= V-)<,-P[C^t], . = l,2,...,iV, 
with bounded 4th order moment. (In fact, ^ 1, for all and i.) Hence, 

oo oo 

E[ j^iY^'^'it) - E[ Y^it) ]r] = Y: n (Yi'^'it) - E[ r^^)(t) ]r ] < oo, 

N=l N=l 

which implies 

Y^!^\t)-E[Y^!^\t)]-^0, a.e., as ^ oo. 
On the other hand, definition of Poisson random measure implies 

-I ^ 1 ^ /-oo 

nY^'^\t)] = ^J2F[rS^^t] = l-J2il-e-^^'^^^'^^^^ / e-Xf\ds), 

which converges to (fTTj) by assumption. □ 

Since by Proposition 11.11 we have almost sure convergence at each time t, we have 
almost sure convergence for all rational number times simultaneously. By definition, 
ycit) and Y^^\u){t), w G fi, are non-decreasing in t. Hence, if yc{t) is continuous, we 
have almost sure convergence as a function in t. 



Corollary 1.2. In addition to the assumptions in Proposition assume that \t 
is continuous in t with respect to the topology of weak convergence. Then for almost all 
sample u E Q, Ylf^\uj) : M+ — [0, 1) defined by ( ITOl) converges pointwise int as N ^ oo 
to a deterministic function yc '■ IR+ — > [0, 1) defined by f|TT]) . O 
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Proposition II . II is a generalization to inhomogeneous case of [121 Proposition 2] for the 
(homogeneous) Poisson process. The correspondence with Aj in Proposition 11.11 and A in 
[I2] is given by At((0, ct]) = A((0, c]). ([9]) imphes that At is the asymptotic distribution of 
the expectation of number of jumps to rank 1 for each particle in the time interval (0,t]. 

Consider a joint empirical distribution /i*^^^ of intensity measure p[^^ and scaled posi- 
tion l^/^'' of the stochastic ranking process: 

1 ^ 

(12) /^!'^^ = ]v5ZV\^r'W)' ^ = °- 

i=l 

p^f^^ , G N, are random variables whose samples are distributions on the product space 
A^(]R+) X [0, 1) of space of Radon measures A^(M+) and an interval [0, 1) C IR+. 

We consider the standard vague topology on A^(M_|_), that is, a sequence {pn} C 
A^(]R+) converges to p G A^(M+) if and only if 

(13) hm / f{s)pr,{ds)= [ f{s)p{ds), 

Jr+ Jr+ 

for all continuous function / with compact support. Since M+ is a Polish space, i.e., 
complete and separable metric space, so is A^(R+) [21 Theorem 31.5], and consequently, 
MiR+) X [0, 1) is also a Polish space [3 Example 26.2]. 
Assume that a sequence of initial configurations 



1 ^ 



(AT) 

- ^-(pI" -1)/N) 
i=l 



converges weakly as A^ — )■ cxd to a probability distribution pq on A^(]R_|_) x [0, 1). Then, 
in particular, 

^^^^ A(^)(rfp):=p(^)(dpx[0,l)) = -^5^(.,(dp) ^ A(rfp):=po(^^px [0,1)), 
weakly, as A^ — )• oo. 

Note also that A^^-* in (|9]) has an expression 

(15) aS^) = / 5,((o,])AW(rfp). 

Jm{r+) 

We shall generalize f lT^ and define, for ^ s ^ t. 



(16) Xif = / 5p((.,t])A(^)(^P)- 

Jm{r+) 
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Theorem 1.3. Assume that /Xq^^ — )■ /iq weakly as N ^ oo for a probability distribution 
/io on A^(]R_|.) X [0, 1). Assume that for each {s,t) satisfying t^s^O, 



(17) A^^^ Xs,t ■= / (5p((s,t])A((ip), weakly as N ^ oo, 

Jm{r+) 



where A is as in (fT4|) . Then for any t > 0, and for almost all sample u & Q, the distribu- 
tion fi[^\u)) converges weakly to a non-random probability distribution fit on A^(]R+) x 
[0,1). 

fit has a following expression in terms of U{dp,y,t) := fit{dp x [y,l)). 

(18) U{dp,y,t):=fit{dpx[y,l)) = ^ -p{{o.t]) rrf^ V +^ m< /i 

[ e ^i>U [dp, y{y, t), 0) yc{t) ^ 2/ < 1. 

Here, tQ{y,t) is the inverse function with respect to to of 

(19) yAito,t) = l- [ e-^«*-*'"*l)A(dp), O^to^t, 
namely, 

(20) to(y, t) = mi{s G [0, t] ; 2/^(s, t) ^ y}, 
and y{y,t) is the inverse function with respect to y of 

(21) yB{y,t) = l- [ e-^«°'*lVo(rfpx [2/,l)), t ^ 0, ^ y < 1, 

Jm{r+) 

namely, 

(22) y{y, t) = inf {x E [0, 1) ; ysix, t) ^ y}. 

O 

Note that yc(t) = yA(t, t) = yB{0,t). Note also that, as will be evident from the proof 
of Theorem 11.31 in Section [2] for ^ y ^ yc{t), the assumption p[l^^ — )■ po can be replaced 
by a weaker assumption Aq^^ A for ^ y ^ yc{t) . 

In contrast to Proposition II. H we do not have a result analogous to Corollary 11.21 for 
Theorem II. 3^ because we can expect no monotonicity for p['^^ . If we impose additional 
conditions, we may go further and prove almost sure convergence as sequences of processes 
on a finite time interval [0,T], both for Y^f^^ — )■ yc and p^^^ — )■ p. See Section H] for 
statement (Theorem 14.11) and proof. 

The structure of the explicit limit formula (ITSIl . in particular, the appearance of the 
inverse functions to of Va and y of ys , can mathematically be understood through a 
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system of partial differential equations, which is a generalization of that in [13] . To avoid 
notational complication, consider the case that the limit distribution A is supported on a 
discrete set: A = Y.^ra,5p^ . Then ([18]) implies, for Ua{y,t) := fJ^t{{Pa} x [y, 1)), 



(23) t) 



f/„(y(y,t),0)e-''«°'*l) yc{t)^y<l, 



where to and y are inverse functions, respectively, of 

(24) l/^(to,t) = l-5^r,e-^'^((*"*«'*l), 

a 

and 

(25) yB{y,t) = l-J2Ua{y,0)e-'-^^'^'^\ 

a 

defined by (|20]) and (|22]). 

Theorem 1.4. Let k be a positive integer, and for each a = 1,2, . . . , k, let he a 
positive constant, Wa '■ IR+ ^+ a measurable function satisfying Wa{t) > 0, t ^ 0, and 
Ua '■ [0, 1) — i- M+ a non-negative smooth strictly decreasing function, satisfying 

k k k 

(26) ^r^ = l, ^r pit) <oo, t^Q, and ^up{y) = I - y, Q y < I. 

/3=1 /3=1 p=l 

Then an initial value problem for a system of partial differential equations 

^27^ ^ ^) + E ^) ^ ^) = (^) ^) ' 

{y,t) e [0,1) xM+, a = l,2,...,A;, 
wrf/i a boundary condition 

(28) f/„(0, t) = r^, t ^ 0, a = 1, 2, . . . , A;, 
and initial data 

(29) ?7„(-,0) =u«, a = l,2,...,fc, 

/ias a unique time global classical solution, whose formula is given by (123]) with 

(30) pa{{s,t\) = / Wa{u)du and Ua{y , 0) = Ua{y) ■ 
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As in [131 §2], ( pTj) is solved by a method of characteristic curves, and Vb, and yc 
turn out to be the characteristic curves for ([27]), which mathematically explains how the 
inverse functions of these functions appear in the solutions. 

For the homogeneous case {Pi^\{0,t]) = w^^h), Theorem 11.31 reduces to [121 Theorem 
5] (with slightly weaker assumption on fiQ, A, and Xt, and with stronger convergence in 
(fi, J-", P), thanks to technical refinement in the proof), and Theorem 11.41 reduces to [131 
Theorem 1]. Motivation for extending the previous results to the present case arises both 
from mathematical and application point of view. 

Mathematical: The model is a natural extension of [12] , with (homogeneous) Poisson 
processes in the formulation of |12] generalized to (inhomogeneous) Poisson ran- 
dom measures in (|2]) or ([5]). Also, as seen from Theorem 11.41 the system of PDE 
corresponding to the limit distribution is a natural extension of that considered in 
[13], with constant coefficients Wa in [13] generalized to time dependent coefficients 
Wait) in (pTj) . On the other hand, the space on which /i^ is defined becomes large; 
/it considered in \12l is a distribution on x [0,1), whereas fit in Theorem 11.31 
is on A^(M+) X [0, 1). Hence it is necessary to extend the definition of the model, 
compared to [121 113] ■ 

Application: The model has successfully been applied to statistical explanation of rank- 
ing data at an online bookstore Amazon.co.jp [131 US] and data of list of subject 
titles at a collected bulletin board 2ch.net [13]. These data arise as results of social 
activities, hence it is inevitable that the data have day-night difference in their time 
dependence. This motivates considering the inhomogeneous cases from an applica- 
tion side. 

Note that we directly see from (|2]), the Markov property 



where we put (A) = \A + u). In practical application, this property enables us 
to shift the time origin t = to the time that a particle we observe jumps to the top, 
namely, we may set x\^\q) = x[^^ = 1, by adjusting the 'clock' for the intensity mea- 
sure accordingly. This motivates our formulating the model in terms of Poisson random 
measures, even though in Proposition 11.11 we apparently do not use Markov properties. 

Note also that if xf^^ = 1, then up to the first jump of i to the top, namely, for 
t < t\^\ comparison of (|5]) and (JH]) leads to 




xf)(t) = xr 



W + 1 
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because, if x\ =1 then, x\, > x\ for all i' ^ i. Therefore, in practical application, we 
may proceed with observing a trajectory (time development) of a single particle, putting 
the time of its first jump to top as t = and observing until its next jump to top, and 
then apply Proposition ll.il or Corollary 11.21 [13| 

The plan of the paper and a brief description of the role of the authors are as follows. 
In Section [2] we prove Theorem 11.31 and we prove Theorem 11.41 in Section [31 In Section HI 
we state and prove Theorem 14.11 time-uniform results corresponding to Proposition 11.11 
and Theorem 11.31 The core structure of the present work, including basic properties of 
the stochastic ranking process which are essential for the proofs of these results, are based 
on collaboration of K. Hattori and T. Hattori. In extending the previous results for the 
convergence of empirical distribution on R_|_ x [0, 1) to A^(]R_|_) x [0, 1), where M.(R^) is a 
space of Borel measures, we have to reformulate the process using Poisson random mea- 
sures and provide abstract measure theory result Lemma [2?H for which collaboration with 
Hariya is crucial. Convergence result as measure valued processes developed in Section H] 
is achieved by collaboration with Nagahata. Also, various technical refinements, implying 
in particular stronger convergence with less assumptions for the uniform intensity case 
[T2] , are results of the collaboration of these 4 authors. In Section [5] we consider a simple 
case where the intensities of the Poisson random measures have a common time depen- 
dence, and prove another scaling limit for the particle trajectory, corresponding to a time 
change with respect to the intensity. This is a result of collaboration of T. Hattori, Hariya, 
Kobayashi, and Takeshima at Tohoku University, and provides a mathematical result of 
scaling limit with time changes, as well as a practically useful formula in applying the 
present results to online rankings. A practical method based on this mathematical result 
is partly checked by actual data obtained at 2ch.net in the master theses of Kobayashi 
and Takeshima (unpublished). In Appendix, we give remarks to be kept in mind when 
applying our results to practical data through statistical analysis. 
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2 Proof of Theorem I1.3L 



Throughout this section, we assume that the assumptions of Theorem 11.31 hold. 
We first note the following rather technical generality. 

Lemma 2.1. Lett > 0. //, for eachy G [0, 1) and for each bounded continuous function 
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g : A^(]R+) — !■ M, there exists Cl with P[Cl] = 1 such that 



1=1 

(31) f /" (7(p)e-''«*-*°(^-*)-*])A(rfp) O^y^ycW, 

-) 

^7(p)e-''«°'*lVo(rfpx [y(y,t),l)) yc(t) < 1, 

IM(M.+ ) 

holds on Cl, then the claim of Theorem \l.'d\ holds for this t. O 

The point here is that Q may depend on y and g, while Theorem 11.31 claims the 
existence of a sample set, independently of test functions. 

We make use of the results in [2| Exercises 30.3, 31.2] for a proof of Lemma [2.11 Note 
that A^(M+) is not locally compact, while local compactness is assumed in the relevant 
results of the reference. We prepare the next Lemma to fill the gap. 

Lemma 2.2. There exists a countable subset T = {/„ ; n G N} of uniformly continuous 
functions fn '■ A^(]R_|_) x [0, 1) — t- R, such that if for each fn^T 

(32) lim / fnip,y)i^N{dp X dy) = fn{p,y) u{dp x dy) 

holds for a sequence of Borel probability measures z/^v and a Borel probability measure v 
on A^(R+) X [0, 1), then z/jv — ?• v, weakly as N ^ oo. O 

Proof. We noted below ( IT3|) that A^(M+) x [0, 1) is a Pohsh space. Note also that 
there exists a coutable set of continuous functions {e„ : M+ — )■ M ; n G N} of compact 
support, such that 
(33) 

c?((pi,?/i), (P2,y2)) = bi - I/2I + ^2""(1 A / en{s)pi{ds)- en{s) p2{ds 

(p„y,) e A1(M+) X [0,1), z = l,2, 

defines a metric d compatible with the topology we are considering [21 (31.4)]. 

Denote a set of sequences by M°° = {x = (xi, X2, . . .)}, and define a metric d' on 
M~ X [0, 1) by 

(34) d'iixi, yi), {X2, 2/2)) = \yi-y2\ + J2 ^ 1^1-" - ^2,n|) 

nGN 

where Xi = {xi^i, Xj,2, . . ■), i = 1,2. We have a natural one-to-one map l = {li,L2, ■ ■ ■ , lq) : 
M{R+) X [0, 1) ^ X [0, 1) defined by 



(35) L{p,y)n= I Cnis) p{ds), n eN, and i{p,y)o = y. 
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Put 



(36) = X [0,1)) C M~ X [0, 1). 

Then ( l33ll . (jSlj) and ( l35ll imply that t : A^(M+) x [0, 1) — > £" is a one-to-one onto isometric 
map. Since A^(M+) x [0, 1) is complete, E' is a closed set in M°° x [0, 1). 

Let F C A^(R_|_) X [0, 1) be a closed set. Since l is isometric, l{F) is a closed subset of 
E' , and since £" is a closed set in R°° x [0, 1), 6(F) is a closed set in R°° x [0, 1). Hence, 
if a sequence of probability measures vn ° ""^ on ]R°° x [0, 1) converges weakly as — i- oo 
to u o L~^, then 

lEK i^^(F) = Ih^ z/^ o i"^(i(F)) ^ o l-\l{F)) = v{F), 

which implies un — )■ z/, weakly as — t- oo. Thus the conclusion of Lemma [2.21 is reduced 
to a weak convergence un o u o l^^ on x [0, 1). 

For each /c G N define a projection to finite dimensional space tt^ : M.°° x [0, 1) — )■ 

X [0, 1) by 

(37) Hk{x) = (xi, X2, . . . , Xfc, y), X = (xi, X2, . . . , y) G M°° X [0, 1). 

Then z/^v o 6^"*^ o and z/ o o jif^^ are probability measures on M'^' x [0, 1). Note 
that a Borel probability measure on Polish space is a Radon measure [21 Theorem 26.3], 
and that the vague convergence of probability measures to a probability measure on M.^ 
is equivalent to the weak convergence P, Theorem 30.8]. Since M.^ x [0,1) is a locally 
compact Polish space, there exists a countable subset Tfc = {fk,i ] i E N} of continuous 
functions fk,i ■ R*^ x [0, 1) — )■ R with compact support, such that if for each fk^i G Tk 



(38) 



lim / fkA^)^N°L ^oTx^^{dz)= I fk,i{^)i^oi ^o7r^^(dz) 

J]R*x[0,l) JR'=x[0,l) 



holds, then o l ^oti^^^uol ^ott^^, weakly as — ?■ oo [2, Exercises 30.3, 31.2]. 
Let 

T = [j{fk,i oTXkOL: A^(R+) X [0, 1) ^ R ; fk,i G Tfc}, 

k&i 

be the T in the assumption of Lemma 1221 Since fk,i, vr, l are continuous, the functions 
in T are continuous. Note further that since fk,i is of bounded support, the functions in 
T are uniformly continuous. Since a countable union of countable sets is countable, T so 
defined is a countable set. With this choice of T, the assumption ( l32l) . with a change in 
integration variable 2; = vr^ o t{p,y), implies 



lim / fk,i{^)'^N°i^ ''°'^k\d.^)= 1™ / 



fk,i o TTfc o ^(P, y) i^Nidp X dy) 

c[0,l) "^'^ JM(R+)x[0,1) 

/ o TTfc o t(p, y) u{dp X dy) = / /fc,i(2;) z/ o o vr^ ^(d^), 

'X{M+)x[0,l) jR'^xIO.l) 



11 



for all k and i, which, as noted below (1381) . implies o t^^ o tt^^ — > z/ o l^^ o tt^^, weakly 
as N ^ oo, for all k. This implies that as measures on R°° x [0, 1), i^n ° i-^^ v o i^^, 
weakly as — )■ oo [3l §2 Example 2.4]. As noted in the paragraph between ( 1361) and ( |37j) . 
this further implies vn — > ^, weakly as — t- oo. □ 

Remark. We could alternatively make use of separability of A^(]R+) directly to obtain 
a countable set T, following the discussion in [201 §1, Remark 4.17, and remark after 
Corollary 9.3]. O 

Proof of Lemma 12.11 Let T be as in Lemma 12.21 If there exists, for each n e N, 
Cln (Z VL such that fl32|) holds for G f2„ and P[r2„ ] = 1 holds, then VL' := f^^^^ f2„ 
satisfies P[ ] = 1 and (l32l) holds for all u & Vt' and /„ G T, which, with Lemma [2.2[ 
implies Theorem 11.31 

Let d be the metric on A^(M+) x [0, 1) as in the proof of Lemma I2l2l Let /„ G T. Since 
fn is uniformly continuous, for any e > there exists 5 > such that for any pi,p2 G 
A^(M+) and yi,?/2 e [0, 1), d{{pi,yi), (p2,Z/2)) < 5 implies |/„(pi,yi) - /„(p2,l/2)| < e. Let 
be a positive integer greater than 1/6 and put 

k-l 

(39) fn,k{p,y) = ^fn{pJ/k)X[l/k,{l+l)/k){y), 

1=0 

where Xla,b){y) = lifa^y<b and otherwise. Then for each p G A^(M+) we have 

sup \fnip,y) - fn,kiP^y)\ < 
2/G[0,l) 

Therefore, limfc_^oo fn,k = /n uniformly on A^(R_|_) x [0, 1). Noting that 

Xll/k,{l+l)/k} = X[l/k,l) - X[(Z+l)/fc,l) , 

we see from (j39ll that has an expression 



fn,k{P^y) = ^9n,k,lip)X[l/k,l)iy), 



=0 



where gn,k,i '■ A^(K+) — )► M is bounded continuous. 

Therefore, if (!3T|) holds, then using the definition (fT2|) and the explicit formula (fTS!) 
claimed in Theorem 11.31 we see that there exists Qn,k satisfying P[ Qn,k ] = 1 and 

lim / fn,k{,p,y) P?\dp X dy){u) = / fn,k{p,y) Pt{dp x dy) 

J M{M.+ )x[0,l) ^A^(M+)x[0,l) 

ifl G (ln,k ■ Hence, f2„ = HfeLi satisfies P[ ] = 1 and (15^ holds for w G ll„ . □ 
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In view of Lemma 12. we fix {y,t) and a bounded continuous function g, in the 
remainder of this section. Since g is bounded, there exists a constant M > such that 

(40) \gip)\^M, peMiR+). 

Since the jump times {tI;^^} are independent. Proposition 11.11 is proved in a straightfor- 
ward way. In contrast, {y}'^^} appearing in the left hand side of ( !3TI) are dependent, and 
moreover, the non-linearity in (127|) indicates that the dependence cannot be neglected in 
the limit oo. A strategy, inherited from the proof in [12], is to (i) choose a nice 

quantity defined as a sum of independent random variables in such a way that the quan- 
tity converges to the right hand side of (13T!) . and (ii) show that the difference between 
the chosen quantity and the left hand side of (13T1) can be shown to disappear in the limit, 
using the properties of the model. We state these two steps explicitly in the following two 
Lemmas, respectively. 

Lemma 2.3. The following hold. 

(i) ForO^y^ycit), 

1 ^ r 
(«) ^ E l.'«((.-..M.l)>o ^ / 9(p)(l-e-«'-'"'«*'»)A(rfrt, 



almost surely as N ^ oo. 

(ii) Foryc{t) ^y<l, 
(42) 

N 

E 9ipf^) ^ / 9{P) e-^«°'*l) poidp X [y{y, t), 1)) 

almost surely as N ^ oo. 
Lemma 2.4. The following hold. 

(i) ForO^y^ycit), 

1 ^ 

i=l 

almost surely as N ^ oo. 

(ii) Foryc{t) ^y<l, 

1 ^ 

i=l 

almost surely as N oo. 
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Proof of (j3l]) assuming Lemma |23] and Lemma [231 For the case ^ y ^ yc{t), 
OD, (H, dUD, and (gSD imply 

TV 



■p{(t-tofe,t),t]) 



A(t;p) 



i=l 
1 ^ 

^ J^^^l/'i'^ ) (- lyW(t)<, + l.f )((i-io(y,t),t])>o^ 
■t=i 

/ 1 ^ TV 

^ i=i 



^(p)(l-e-''«*-*»(J''*)'*l))A(rfp) 



A4{R+) 



TV 



= ^^I^|lv('V),.,_-l.,(iV) 



+ 



+ 



1=1 

^?(p)AW(rfp)- / ^?(p)A(rfp) 

M{R+) Jm{R+) 
N 



AT^-- ^ ^'^^^ ^((t-to(.,t),t])>o ^^^^^^ 



g{p) (1 -e-''«*-*«(f'*)'*l))A(dp) 



— 0, a.s., as — > oo, 

which proves f l3T|) for ^ ?/ ^ Vcit)- 

Similarly, for the case yc{t) < I, (HQ]), (I42]), and (HI) imply 



j=l 



(7(p)e-''«°-*])po(rfpx[y(y,t),l)) 



A1(K+) 
1. 



g(p)e-'«''''n^,„{dpxly{y,t),l)) 



1=1 



^(x(^)-l)/TV^j/(s,,t),r(f >t 



< 



1 ^ 

mIv |i 



TV 



^ 5Z^7(pf ^) l(.f)-l)/TV^^(,,t),.(f >t 

^?(p)e-^«°'*Vo(rfpx[y(y,t),l)) 



which proves ( ETj) for yc(t) ^ y < 1. 



— 7- 0, a.s., as — )■ oo. 



□ 
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Before proving Lemma 12.31 and Lemma 12.41 we prepare a couple of random variables 
which converge as — t- oo to ?/a in ( IT9|) and ub in (EI]) . The following Lemma 12.51 is used 
in the proof of Lemma 12.41 and the proof of Lemma 12.31 is similar to that of Lemma 12.51 

Lemma 2.5. (i) For ^to define 



1 ^ 

(45) yrito,t) = ^j:i 



N A^^yf\{t-toA)>o■ 



Then Y^\tQ,t) yA{to,t), almost surely as N ^ oo. 
(ii) For t ^ and ^ ?/o < 1 define 

(46) >^r(yo,t) = yo + ^ E K^-h 



i; {x\^^-l)/N^yo 

Then Yl^^\yo,t) — > yB{yo,t), almost surely as N ^ oo. <> 



Proof. As in the proof of Proposition II. ![ a strong law of large numbers implies, 
almost surely as — t- oo 

Yf\to,t) - E[ Yf\to,t) ] ^ and ^(yct) - E[ Yj,''\yo,t) ] ^ 0. 

On the other hand, (fT6!) and (fT7|) imply 

1 ^ 

hm E[ Yf\to,t) ] = lim - ^^(1 - P[ ^((t - t,,t]) = ]) 

i=l 

, N „ 

= lim _y(l-e-''''''«*-*'"*])) = l- lim / e-''«*-*«'*])A(^)(rfp) 

POO p 

1 - hm / e-^xl^ids) = 1-1 e-^«*-*-*l)A(dp) = y^(to, t). 



Similarly, 



M{R+) 



lim E[ ^2/0, t) ] = I/O + lim J] P[ ^ t 

j; {a;<'^'-l)/Ar^yo 

= ^o+liml Y: (l-e-^'"«°'*^)) 

N-¥oo iV 



1/0 + 1 - yo - lim / e-^«0'*Vi'^Hrfp x [yo, 1)) 
e-''«°'*lVo(^^px [1/0,1)) =l/B(z/o,t). 



□ 
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Proof of Lemma 12.31 The proof is a repetition of the proof of Lemma [23| by replac- 
for (il. ' ' □ 



The proof of (13T!) now will be complete if we prove Lemma 12.41 which is proved in a 
similar way as the corresponding part in [12]. 

Proof of (|43]) for ^ Vcit)- Note that VAitoit) of (US]) is non-decreasing in to 
and t, with ?/^(0,t) = and yAit,t) = yc(t), and by assumption of the Theorem 11.31 is 
continuous. Hence 

(47) VAitoiy, t),t) = y, O^y^ yc{t), t ^ 0. 

Lemma 12751 therefore implies that there exists Qa C Q, satisfying P[fiyi ] = 1, such 
that 

(48) hm Yi''\to{y,t),t){u) = y, ueQA- 

Fix u G Qa arbitrarily. The definition of the stochastic ranking process and ( l45ll im- 
ply that i'l'^\{t — to{y,t),t]){u) > 0, if and only if Y}'^\t){u) is on the top side of 
Yi''\to{y,t),t){ooy, Y}''\t){u;) < Yi''\to{y,t),t){oo). Therefore 



1 ^ 



(49) 



i=l 

Note that the definition of Y}^\t) in ([6]) implies that it takes values in {k/N ; k = 
0, 1, . . . , - 1}. Hence (gH]) implies 

I ^ 

^ hm i- X {N\Yf\to{y,t),t)iu)-y\ + 1) = 0. 

Af— i-oo iV 

The relations dH]) and (jSO]) imply (gS]). □ 



Proof of f|44]) for ycit) ^ y < 1- ysiy^t) of fl2T]) is non-decreasing in y and t, with 
ysi^^t) = yc(^) and yB{l—,t) = 1 — 0, and by assumption of the Theorem [L3l is contin- 
uous. Hence 

(51) yB{y{y, t),t) = y, ydt) ^ y < 1, t ^ 0. 
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Lemma [2.51 therefore implies that there exists fi^ C f2, satisfying P[ fi^ ] = 1, such that 
(52) lim Yi,'^\y{y,t),t){uj) = y, oo e . 

iV— !-oo 

Fix u E Qb arbitrarily. The definition of the stochastic ranking process and P6l) imply 
that {x\^^ — 1)/N ^ y{y,t) and t^'t^\uj) > t hold together, if and only if Y}^\t){uj) is on 
the tail side of Y^^\y{y,t),t){u); Yl^\t){u) ^ Y^^\to{y,t),t){u). Therefore 



1 ^ 



(53) 



As in the proof of (143|) . F/^''(t) takes values in {k/N ; A; = 0, 1, . . . , A^ — 1}, which implies, 
with dSlD, 

1 ^ 

^ lim X (Ar|Fi^)(y(j/,t),t)(a;) - y| + 1) = 0. 

N-^oo iV 

The relations ([Ml) and (|54]) imply (gl]). □ 
This completes the proof of Lemma 12.41 hence of Theorem II. 3[ 



3 Proof of Theorem II. 4L 



To prove Theorem II. 4[ we apply a standard method of characteristic curves. 

First, assume ^y ^ Ucit) = yAit,t). Let ti ^ 0, and consider an ordinary differen- 
tial equation for a characteristic curve intersecting (0,ti), defined by 

.... -^(t) = ^w^(t)?7^(y(t),t), a = l,2,...,k, t^t,, 

(55) at 

yih) = 0. 

Put 

(56) ip^it) = Ua{y{t),t), a = l,2,...,k, t^h. 
Then ([56D, ([27D, and ([55]) imply 

(57) ^(t) = -w^it)UMt),t) = ~w^{t)^^{t), 
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which, with y(ti) = in ( 155|) . has a unique solution 



(58) y^a{t) = Ua{0,ti) exp(— / Wa{u)du) = Ta exp(— / Wa{u)du), 

Jtl Jtl 
where we also used (128|) . Substituting (l56l) and (158!) in (155|) . we have 



dy 



dt 

which, with y{ti) = 0, has a unique solution 



(t) = ^ w,3{t) rp exp(- / wp{u)du) 



(59) ?/(t) = J^r;3 (1 - exp(- / wp{u)du)) = yAit - ti,t). 

where we also used X]^=i = 1 in (1261) and (1241) with (1301) . in the last equality. The 
assumptions for Wa in Theorem 11.41 imply that yA{to,t) is strictly increasing and differ- 
entiable in to, satisfying ?/^(0,t) = and yAit,t) = ycif)- Hence there exists a unique, 
strictly increasing, differentiable inverse function to = to{y,t), taking values in [0,t], sat- 
isfying 

yA{to{y,t),t) = y, O^y^ycit), t^O. 
This, with ([56]), dSH]), and ([59]), implies 

Ua{y,t) = exp{- Wa{u)du), 

Jt-to(y,t) 

which proves ( l23i) for ^ ?/ ^ ycif)- 

Next, assume ycif) = yB{0,t) ^ y < 1. Let ^ yo < 1, and consider an ordinary 
differential equation for a characteristic curve intersecting {yo,0), defined by 

d ^ 

^gQ^ ^(t) = ^w;;3(t)f//3(l/(t),t), « = l,2,...,fc, t^O, 

1/(0) = yo ■ 

Put 

(61) (^^(t) = f/«(?/(t),t), a = l,2,...,fc, t ^ti. 
Then (|6T1) . fl27]) . and fl60l) imply, exactly as for the case ^ ycit), 

(62) ^(t) = -W7,(t)^,(t), 
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which, with y{0) = yo, has a unique solution 



(63) (fait) = Uaivo) exp(- / Wa{u)du), 

Jo 

where we also used ( l29l) . Substituting ( l6Ti) and ( l63i) in ( l60i) . we have another differential 
equation for y{t), which, with y{0) = yo, has a unique solution 

(64) y{t)=yBiyo,t), 

where we used ^^=iM^(y) = 1 — y in ( l26l) and ( 125|) with ( l30i) . The assumptions for Mq, 
in Theorem II .41 imply that |/b(?/, t) is strictly increasing and differentiable in y , satisfying 
yB{0,t) = yc(t) and yB{l—,t) = 1—. Hence there exists a unique, strictly increasing, 
differentiable inverse function y{y,t), taking values in [0, 1), satisfying 

yB{y{y,t),t) = y, yc{t)^y<l, t^O. 

As in the proof for y ^ ycif)i this, with (!6T]) . (163|) . and f l64|) . implies (123|) for yc'('^) = V < 1- 
This completes a proof of Theorem 11.41 □ 

4 Scaling limit results uniform in time. 

Let T > and 

(65) X = {rj^^ : [0,T] ^ R+ ; i = 1,2, . . . , N, N e N} 

be a set of continuous functions on [0,T] defined by r^^\t) = p\^\{0,t]), t ^ 0. Note 
that since we assumed in the beginning that p^-^"* is continuous, r^^^ is continuous. In this 
section, we prove the following. 



Theorem 4.1. LetT > 0. In addition to the assumptions in Proposition \l.l\ assume 
that a set of continuous functions I defined by ( l65l) is uniformly equicontinuous ; namely, 



(66) limsup sup \r{s) — r(t)\ = 0. 

SiO s,te[0,T]; \s-t\S5 

Then, Y^^^ of ( fTOj) converges almost surely to yc of ( fTTl) as N ^ oo, as a sequence in the 
space of continuous functions on [0, T] with supremum norm: 



(67) 



N 



lim sup \Y^''\t)~yc{t)\ = 



>oo 



i6[0,T] 



Assume next that all the assumptions of Theorem \1.3\ and (l66l) hold. Assume also that 
a set of functions 

J = {r: [0,T] ^ M+ ; r{t) = p{{0,t]), t G [0,T], p G suppt A} 
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is uniformly equicontinuous, and that for i/a of (fT9|) and ys of (pT]) . yAit — ti,t) and 
ysiy^t) are equicontinuous in (ti,t) and {y,t), respectively. Then, jj,^^^ of (fT2!) converges 
almost surely to fi. of (|T8l) as N ^ oo, as a sequence in the space of probability measure 
valued functions fi. : t fit with supremum norm. O 

Proof. First we assume that the assumptions of Proposition II. II and ( !66l) hold. Note 
that dH) imphes that, for i = 1, 2, . . . , A^, = 1, 2, . . ., 

-rmt])-rf\t) = .rmt])~pf\m), t ^ o, 

is a martingale up to fixed time T. Note also that ([3]) implies 

(68) ^f^((0,tAT;(f]) = 
Hence 

(69) lyf Ht) := V.)<,-rf )(t Ar(f ), t G [0,T], 
is a bounded martingale. This with ffTOl) further implies that 

TV Af AT 



is also a bounded martingale. Using Doob's inequality, independence of {t^^ 
1, 2, . . . , N}, and iW^i^^T)! ^ 1, we have 



(iV) 



E 



0<t<T V^V 



i=l 



44 

- 34 



< 



With an argument similar to that in the proof of Proposition 11.11 



(70) 



sup 

0<i<T 



1=1 



0, a.e., as N ^ oc. 



On the other hand, for each ^ t ^ T, as in the proof of Proposition [TTTl independence 
and boundedness of rl^\t A t^'t^^), i = 1,2, . . . , N , imply 



N N 

(71) ^j:rr\t^r^P)-^j:nrritAr^P)]^0, a.e.,as N ^ 00, 



i=l 



i=l 



and 



E[ rf\t A T^P) ] = E[ ] - E[ Wr'it) ] = E[ 
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implies 



(72) 



1 ^ 

Af— )-oo iV 



i=l 



lim E[r^^)(t)]=yc(t). 

W— )-oo 



Since Y^^\t) is non- decreasing in t, and yc{t) is its pointwise limit, it is also non- 
decreasing. As in the case of Corollary \1.2\ (I7T1) and fl72|) imply that, with probability 
one, 



(73) 



1 ^ 
1=1 



Since {'"j-^'*} is equicontinuous, ( 17S|) implies that yc" is continuous on rationals, and the 
monotonicity of yc proves that it is continuous on [0,T]. 

By assumption of equicontinuity and the convergence (I73|l on a dense subset of [0,T], 
it follows that the convergence is uniform: 



(74) 



sup 

t€[0,T] 



i=l 



— )■ 0, a.e., as iV — )■ oo. 



The equations ([TO]), (EH]), (|70]), (l74|) prove (|67|). 

In the remainder of this section, we assume that the assumptions of Theorem 14. ll hold. 
To prove uniform convergence of jj,^^^ , we first prepare t- uniform version of Lemma 12.11 

Lemma 4.2. //, for each y G [0, 1) and for each bounded continuous function g : 
A^(M+) — M, there exists Q with P[Q] = 1 such that, for each u ^ Q, 



g{p) p!f\dp X [y, l)){uj) - I g{p) pt{dp x [y, 1)) 

M{R+) Jm{R+) 



(75) lim sup 

then p[^^ converges to pt uniformly in t G [0,T] as N ^ oo, almost surely. 



0, 



O 



Proof. Let T = {/„ ; n G N} be as in the proof of Lemma [2.11 and for probability 
measures p and z/ on A^(]R+) x [0, 1), put 



7r(/i, z/) 

oo 



n=l 



fnip, y) pidp X dy)- / /„(p, y) iy{dp x dy) 

X(R+)x[0,l) J M{M.+ )x[0,l) 



All 



Then vr is a metric on the space of probability measures on A^(]R+) x [0,1), and the 
convergence with respect to vr is equivalent to convergence (132|) for each /„ G T. Hence, 
as noted just below (1321) . it is equivalent to weak convergence of the probability measures 
on M{R+) X [0,1). 
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Now assume that (!75|) holds. Then following the arguments of the proof of Lemma [2?T| 
replacing (I3T]) by (!75|) . we see that there exists Q' G Q such that P[ ] = 1 and 



lim sup 



fn{p,y)fJ'i^\dp X dy){uj) - / fn{p,y)fJ't{dp x dy) 

X(R+)x[0,l) J M{M.+ )x[0,l) 



0, 



for all n G N and u E Q'. Therefore, 



lim sup 7r{fif^\uj), Ht) = 0, a; G Q' , 

which, by the equivalence of convergence in tt and the convergence in the weak topology 
of the space of probability measures on A^(]R+) x [0, 1), implies the almost sure uniform 
convergence in t G [0,T], of to fit- □ 

In view of Lemma [4.21 we fix ?/ G [0, 1) and a bounded continuous function g, in the 
remainder of this section. Note that ( l40l) holds. The assumption A*^^) — )■ A in ( |T4l) further 
implies that for any K > there exists a positive integer Nq such that, for N > Nq , 
(76) 

M 



(P)A(^P) 



^7(p)AW(rfp)- / ^?(p)A(rfp) 



The following Lemma corresponds to Lemma [2.51 

Lemma 4.3. For eachti G [0,T], rj^^ o/dlS]) andyA of satisfy 



(77) 



sup \Y^^\t — ti,t) — yA{t — ti,t)\ ^ 0, a.e., as N ^ oo, 

ie[ii,r] 

r{N) 



d for each y^ G [0, 1), o/ (H^ and ys of satisfy 



(7J 



sup |Y^^^(?/o,i) - 2/i?(yo,i)| 0, a.e., as N ^ oo. 

te[o,T] 



< 



K 



Proof. Define, for z = 1, 2, . . . , A^, = 1, 2 



i.ki ' 



where ki := inf{j ; t^^^ > ^i}- Then just as in the proof of flB7|) . we see that 

vrit) ■.=^.,^,-pT\{t„tAfr\te[t,,T] 

and, with (H5l) . accordingly, 

TV TV 

(79) rf (t - t„t) - j^Y.P^\^'^^'^ = ^ ' ^ [^-^] 



i=l 



i=l 
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are bounded martingales, and we have 



(80) 



sup 



i=l 



^0, tGQn[0,T], a.e., as ^ oo. 



On the other hand, we have with probabihty one. 



1 ^ 

(81) -5^pfH(^l,^Af^('^)]) = ^/^(^-^l,^), tGQn[o,T], as n^oo. 
1=1 

By assumptions of equi continuity and the convergence (18T!) on a dense subset of [0, T], 
it follows that the convergence is uniform: 



^2) 



sup 

te[o,T] 



1 ^ 

i=l 



— 0, a.e., as N ^ oo. 



The equations ([79]), (EOD and ([82]) lead to ([77]). 

A proof of (178]) goes in exact correspondence with that of (1671) . if we directly use the 
assumption of continuity of in place of monotonicity of yc ■ □ 



Corollary 4.4. For each ti G [0,T], 



(83) sup 

te[ti,T] 



1 ^ TV /■ 



^7(p)(l-e-''«*-*l))A(rfp) 



^0, 



almost surely as N ^ oo, and for each G [0, 1), 
(84) 



sup 

te[o,T] 



N 



N 



i=l 



+ ) 



^7(p)e-''«°'*Vo(cipx [yo,l)) 



-> 0, 



almost surely as N ^ oo. 

Proof. This is proved as in the proof of Lemma [4. 3 [ if one notes (140]) . 



□ 



Fix a positive integer K arbitrarily. By the assumptions of Theorem 14.11 of uniform 
equicontinuity of J', yA and yB, and noting that fiQ{Ai(M^) x [y, 1)) = 1 — y, there exist 
a positive integer L and sequences = ti^ < ti^i < ■ ■ < ti^i = T and = ?/o,o < yo,i < 
■ ■ ■ < yQ,L = 1 such that 



(i) for j = 0, 1, 2, . . . , M and s G [tij_i , tij+i], 
(85) / 



g-p{(il.j.i]) _ g-p{*.t])| 



A{dp) < 



'7M(R+) 



K 



23 



where, for convenience we put tij = if j ^ 0, and tij = T ii j ^ L, and also for 
j = 0, 1, 2, . . . , M and 2; G [yo,j-i , Voj+i], 



^7(p)e-^«°'*])/io(rfp X [yo,, 1)) - / ^?(p)e-^«°'*])po(rfp x [z, 1)) 

) ^A^(K+) 

(7(p)e-^«°'*lVo(c?P X [yoj A z, V z)) ^ 

) 



^6) 



where, we put yoj = if j ^ 0, and yoj = 1 if j ^ L, 



< 



K 



(ii) the sequences of functions yA,jif) = yA{.ii — tij) V 0, t), j = 0, 1, 2, . . . , L, which is 
decreasing in j, and ysjit) = yB{yo,j, t), j = 0,1,2, . . . , L — 1, which is increasing 
in j, satisfy 

(87) ^ yA,j{t) - yA,o+i{t) < , J = 0, 1, 2, . . . , L - 1, t e [0, T], 



and 



^ yB,,+i{t) - yBjit) < , J = 0, 1, 2, . . . , L - 2, t e [0, T]. 



Lemma [4.31 and Corollary 14.41 imply that there exists Cl^ C Q, satisfying P[ Cl^ ] = 1, 
such that for all u G Cl^ there exists an integer Nq = Nq{u) such that if N > Nq then 



(89) 



\yA'\t - ti,j,t){uj) - yA,,{t)\ <j-^,te [ti,„T], J = 0, 1, . . . ,L, 



(90) 



\Yli^\yo,j,t){io)-yBAt)\ < ^, J = 0,1,...,L, te[0,T] 



K 



(91) 



and 



(92) 



1 ^ /• 



)A(rfp) 



< 



M 



[ti,„T], j = 0,l,...,L, 



iV 



< J=0,l,...,i.-1, tG[0,T]. 



Now, we shall consider the case ycif) = y and the case yc{t) = y separately. First, 
let yc(t) ^ y, and let j = j(t) be the integer such that 



(93) 



yA,j{t) ^y < yA,j-i{t). 
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Note that yc(t) ^ y implies y = yAito{y,t),t) (see (H?])), with which yA,o(t) = yA{t,t) = 
ycif)i yA,L(t) = yA{0,t) = 0, and monotonicity of yA{to,t) with respect to to imply that 
such an integer j = j(t) exists if ycit) = y- Since yA{to,t) is increasing in to , (ESI) also 
implies 



(94) 

Since (|871) implies 



^ y - yA,j{t) ^ yA,j^i{t) - yA,j{t) < 



K 



with ( 189|) and a similar argument as for ( 150|1 . we have 



(95) 



i=l 
1 ^ 

i=l 



< 



'^Yl''\t)<y-\^^\t)<y^,,{t) 



i=l 



< 



M{y - yAM) + M\Yi''\t - h„t)iu) - y^,,(t)| < 



2M 



Note also that, as in the argument for 

Adding up (!76|) , (l95l) , (19T]) and (!85l) , and using (p6|) and triangular inequality, we arrive 



at 



sup 

te[o,T]; yc{t)^y 



A (dp) 



< 



(97) 



sup 

N 



^(p)(l-e-''((*-*«(^'*)'*]))A(rfp) 



M{R+) 



< 



5M 



for u eQk and > Ao(i^)- 

Next, let ycit) ^ and let j = j(t) be the integer such that 



(9J 



yB,j{t) ^y < yB,j+i{t). 
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With an argument similar as that below (!93|) . such an integer j = j(t) exists if yc(t) ^ y. 
Since ysiyoyt) is increasing in y^ , ( p8|) also implies 



(99) 

Since flHHj) implies 



^ y - ^ yB,j+i{t) - ysjit) < 



K 



with ( 190|) and a similar argument as for ( 154|) . we have 



(100) 



(V 



i=l 



V/-)(t)^,,,(i)-ly(-)(t)^yrfeo.,t) 



< 



M{y - yB,m + M|Fi^)(yo,„t)(c^) - < 



2M 



Note also that, as in the argument for (!53l) . 



(101) 



Adding up fllOOp . fl92|) and fl86|) . and using f llOip and triangular inequality, we arrive 

at 

(102) 



sup 

iG[0,T]; yc(my 

^ 4M 



<?(p)e-^«°'*lVo(rfpx [y(y,t),l)) 



>1(M+) 



for w G fiif and > A'o('^)- 

Combining (p7I) and (11021) . we have 



sup 

te[o,T] 



g{p)fif\dp X [y, l))(u;) 

A^(M+) 

> A'o( 



g{p)pt{dp X [y, 1)) 



< 



5M 



Finally, put Vt = ni^=i ^k- Then P[ ] = 1- Let w G fi. For any e > take an integer 
K such that K > 5M/e . Then u ^ Cl G Cl^ implies 



sup 

te[o,T] 



g{p) pf\dp X [y, l))(tj) - / g{p) pt{dp X [y, 1)) 



A^(IR+) 



^ 5M 



for A^ > A^o(i^)5 which implies ( !75|) . and therefore Lemma [4.21 implies the Theorem. □ 
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5 Case when the intensities have common time de- 
pendence. 

To consider the case where the intensity measure p has a density, denote the set of locally 
integrable functions on R+ by ^1^^,(1]^+). Lj'Qj,(M4.) is a complete separable metric space. 
Let t be a map l : Lj'^^^(]R+) — )• A^(]R4.) which maps w G Lj'jjj,(R+) to the measure on M+ 
with density w determined by 

(103) L{w){{s,t]) = w{u)du, ^ s <t. 

Proposition 5.1. Assume that w'l^'' e Lj'Q^(M+), i = 1, 2, . . . , A^, = 1, 2, . . and 

for each N, put 

i=l 

If there exists a probability distribution A on Ll^^(K^) such that converges weakly to 
A as N ^ oo, then the sequence of distribution A^^\ N = 1,2, . . ., on the set of intensity 
measures A^(M+) defined by A^^-* = A*^^-* o l^^, with l as in f ll03p . converges weakly as 
N ^ oo to A := A o l^^. Moreover, for all ^ s < t, X^^^^ defined by ( !T6|) converges 
weakly as N ^ oo to Xs^t defined by ([T] 



Proof. Let g : M{R+) — )• M be a bounded continuous function on A^(]R+). Then 
the definitions imply 

^(p)AW(cip)= / g{t{w))m{dw). 

Let {wn} be a sequence converging in Lj'Qj,(]R+) to w, and let / : M+ — t- M be a continuous 
function with compact support: f{u) = 0, u ^ k, for some integer k. Then / is bounded: 
\f{u)\^M,ueR+, for some M. Hence 



f{u)Wn{u)du— I f{u)w{u)du 



k 



k fk 

f{u)wn{u)du— I f{u)w{u)du 
Jo 



\wn{u) — w{u)\du ^ 0, n — 7- cxD. 







This holds for all continuous function / with compact support, hence lim„_j.oo iiwn) = iiw) 
in vague topology, which further implies 

lim g{L{wn)) = gii-iw)). 
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This proves that got: Lj_,^(]R+) — )■ M is a bounded continuous function, hence the 
assumption A*^^) — > A imphes 

hm / g{p) A^''\dp) = f g{p)A{dp). 

TV— i>oo 



1+) JM{R+) 

This holds for any bounded continuous function g, which proves A^^^ — )■ A, weakly as 
oo. 

Let t > s > and put b[w] = w{u) du. In a similar way as above, the definitions 
imply 

><iV = / 6b[u,]A^^\dw) and Xs,t = / 6b[w]Hdw). 

Let h : M+ — )■ M be a bounded continuous function. Then the map 

Ll^(M+) h{b[w]) e M 

is bounded and continuous, hence the assumption A^^^ — )■ A implies 
h{w) Xi't\dw) = [ h{b[w])A^''\dw) 
/ h{b[w])A{dw) = / h{w) Xs,t{dw), N ^ oo, 

hence X[^f^ — )■ Xs,t, weakly as — )■ oo. □ 

Proposition 15.11 implies that the assumption (JT71) in Theorem 11.31 is redundant if the 
intensity measures have densities. 

For the rest of this section, we further assume a common time dependence for all 
w^^'^ in Proposition 15.11 Namely, we assume that there exist d G Lj'j^^(M_|_) and positive 
constants 



wf^>0, t = l,2,...,N, N = l,2,... 



such that the intensity measure of the Poisson random measures v^^^^ in the stochastic 
ranking process ([2]) is given by 

(104) pf\{s,t])=wf^ td{u)du, t = l,2,...,N, N = l,2,.... 



As in the proof of Proposition 15. H we have 

Corollary 5.2. Let d G Lj'jjj,(M+). If there exists a probability distribution X on 
such that 



1 ^ 

(105) A^^^ := — ^ 6jN) A, weakly, as N ^ oo, 



i=l 
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then a sequence of probability distributions A^^\ N = 1,2, . . ., on LJq^(]R+) defined by 

r 1 ^ 



converges weakly to a probability distribution A = f^^ 5wa^{dw), as N ^ oo. 

In particular, Proposition holds with Pi^\{s,t]) = w^^^ d{u) du, and yc{t) of 
( ITTj) is given by 



(106) yc{t) = l- / e--^WA(ciu;), 
where 

(107) A{t) = [ d{u) du. 



O 



The formula fll06p is to be compared with the case of the (homogeneous) Poisson 
process in [121 Proposition 2], where we have 



;i08) ycit) = 1-1 e-^'\{dw). 

+ 



A in fllUSj) is the (infinite particle limit asymptotic) distribution of jump rates, while A in 
the case of common time dependence fll06p is the distribution of relative jump rates. 

To study a time change according to the common intensity measure, let us first make 
a heuristic observation. Suppose we could trace the trajectories of n ^ particles 
3ii32i ■ ■ ■ ,jn- The total number of jumps of the n particles in the time interval (0,t] is 
given by 

n 

(109) 5(^'")(t) = 5^z.f)((0,t]). 

i=l 

If n is large (n ^ 1), we expect as a consequence of the law of large numbers, as in 
Proposition 11.11 

(110) 5(^'")(t) ^ Epr((0,t]) = A{t)Z{N,n), 

i=l 

where we put 

(111) Z{N,n) = J2w^^ 



u 

i=l 
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and also used (fTOD and (fTOT]) . Using ffTTOD in (fT06D . we have 
(112) yc{t)^l- [ e-^^^'^'"H^)/^(^'-)A(rf 



Wj. 

+ 



The approximate formula (I112p suggests that, if we perform a time change t' = S^^'"'\t), 
then modulo scaling constant Z{N,n), we recover a formula (11081) for the homogeneous 
case. 

We can put the heuristic consideration which lead to (11121) in a mathematically precise 
form. For t ^ 0, let 

N 

(113) 5(^)(t) = 5^z.f)((0,t]) 

and denote its right continuous inverse by 

(114) s^^\t) = mf{s ^ ; S^^'^s) > t}. 
Let d G Lj'^j^(]R_|_). For simplicity, assume further that 

(115) fi(t) > 0, t^O. 

Then A{t) of (11071) is strictly increasing, and the inverse function is also continuous. 
Theorem 5.3. Let d e Ll^^{R+), and assume (I115p . Put 

N 

(116) Z{N) = J2 
and assume 

(117) lim Z{N) = oo. 

A''— >-oo 

//, as in Corollary \5.2\ there exists a probability distribution A on M+ such that (11051) 

holds, then for each t ^ 

(118) 

Y^^\s^^\Z{N) t)) -> yciA-\t)) = 1- [ e-'"'\{dw), m probability, as N ^ oo, 

Jr+ 

where Y^"^ is defined in (fTOj) . O 
To prove Theorem 15. 3[ we first provide a rigorous version of (IllOp . 



N 
i=l 
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oo. 



Lemma 5.4. Fort ^ 0, 
(119) 
and 

(120) s'^^\Z{N)t) A-\t), m prohahiUty, as N ^ oo. 



-^^^S'^^^(t) A(t), in probability, as N ^ 



Proof. Since by definition u^^\{0,t]) follows the Poisson distribution with expecta- 
tion pf^\{0,t]), we have 

(121) E[ S^^\t) ] = V[ 5(^)(t) ] = A{t) Z{N), 

where V[ ■ ] denotes variance. For e > 0, fll2ip . flll6p . and Chebyshev's inequahty imply 

P[ |5W(t) - E[ 5(^)(t) ]| > Z{N)e ] ^ {eZ{N))-'Y[ ^^(t) ] = , 



which, with f lll7p . imphes 
1 



{S^^\t) -E[S^^\t)]) ^0, in probabihty, as N ^ oo. 



Z{N) 

This, with f lT2T|) . imphes (fTT9|) . 

Next, noting that S^^\t) is non- decreasing in t, (I114p implies 
(122) 

{uen; .(^)(Z(iV)t)(a;) ^ A-\t) + e} C {co e ; ^^^""^ (^"^ W + ^ M ^ 0- 

The assumption flllSp implies that A is strictly increasing, hence, S = A{A^^{t)+e/2)—t > 
0, and 

Wen; + 



C {to E^l; 



Z{N) 

This and ffTTOj) and ([I22]) imply 



'-^/''\A-\t)+'-)-A[A-\t) + ^) 



> 



6}. 



(123) 



lim P[ s(^)(Z(iV)t) ^ A-i(t) + e ] = 0. 



Similarly, 6' = t - A{A-\t) - e/2) > 0, and 

{ioen; s^^\Z{N)t) ^ A-\t) - e} 
C{cuen; TT-U^^^) ( A'\t) -'-)^t} 



Z{N) 



\^S^'^\A-\t)-^)-A[A^\t)-'^ 



Z{N) 
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lim F\s^^\Z{N)t) ^A~\t)-e 

N^oo 



0. 



which imphes 
(124) 

(Imll and (IT2411 prove (IT20D . 

Proof of Theorem 15.31 By triangular inequahty, we have 

\Yf\s'^^\Z{N)t))~yc{A-\t))\ 

^ |yf )(.W(Z(iV)t)) + \Y^^\A'\t)) - yc{A-\t))\. 



□ 



Corollary 15.21 implies that the second term in the right hand side converges to in prob- 
ability as — J- oo, so it suffices to prove that, for all e > 0, 



hm P[ |F^^)(.W(Z(iV)t)) ^ 



(125) 

holds. 

For 6 > put 

(126) fif^={u;G^]; \s^''\Z{N)t){uj) - A-\t)\ < S}. 
Then fll20p implies 

(127) lim P[r]J^^M = 0. 







The definition ( JTOl) of Y^f^^ implies 

\Y^'^\s(^\ZiN)t))-Y^''\A-\m 



(128) 



TV 



Combining f ll26p and f ll28p . we have 
P[|F^^)(.(^)(Z(iV)t))-F^^)(A-i(t))| ^e, fi,] ^P 



TV 



5^ l-i.f 



e(A-i(t)-5,A-i(t)+<5) 



i=l 



Applying Chebyshev's inequality, we further have 

P[|r^^)(.(^)(Z(Ar)t))-F^^)(A-i(t))| ^e, ns] 
< 



N 



ttY^I 1 (iv 

i=l 
1 ^ 



G(A-i(t)-5,A-i(t)+5) 



A^e 
1 
e 
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This, with f ll27p and the assumption f llOSp . imphes 

lEK P[ |F^^)(.W(Z(iV) t)) - Y^!'\A-\t))\ ^ e ] 

N—^oo 

^ ]I^ P[nf'^"]+m, F[\Y^''\s^''\ZiN)t))-Y^'^\A-\t))\^e, nf^ 

\{dw). 



< 1 

~ e 



1 



7V-i>oo 



^-AiA-Ht)-S)w _ ^-AiA'Ht)+S)w 



This holds for all 5 > 0, hence the bounded convergence theorem and the continuity of 
A{t) imply 

lE^ P[ |y^^)(.W(Z(iV) t)) - Yf \A-\t))\ ^ e ] 

N^oo 

< inf - / fe-^(-^"'W-^)'" - e-^(^"^W+^)-) \{dw) 

< - [ hm fe-^(^"'«-^)'" - e-^(^-^W+^)-) X{dw) = 0. 
This proves (11251) . hence Theorem 15.31 is proved. 



□ 



As an explicit example to Z{N) and A, consider, as in [131 E]) the Zipf's law, which 



IS 



(129) wr 



N 



1/6 



a(-) , « = l,2,...,iV, 



for positive constants a and b. For this choice. 



N 



0—1 



6 > 1, 



(130) Z{N) = J2wf ^ = (1 + o(ArO)) ^ ^ 



aA^ / x'^dx = aN log N b= 1, 

il/iV 

oo ^ 

5Z ^ = a^'^'C(l/&) < 6 < 1 



The corresponding N ^ oo weak limit is the (generalized) Pareto distribution, defined 
by 



(131) 



\{[w,oo)) 



w 



w ^ a, 
w < a. 



With the Pareto distribution (fT3T]) for A, (fTT2|) is (for N = n) 
(132) 



xc{t) = Nyc{t) + N - N 

b 

T(-b 



N-b 
N -Ne 



/ gw(t) Y 



s(^)(t) 



'iVi/bCiv(l/fe)- 



^5(^)(t)/(Afi/''Civ(l/6)) I 
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where Cn{z) = ^i=i^"^- The last hne in f ll32p is obtained by integration by parts from 
the second hne, as in [T3], and is suitable for < 6 < 1. Note that the parameter a in the 
Pareto distribution fllSip disappears in the time changed formula f ll32p . 

A Remarks on practical application. 

In fl3[ [H] , the mathematical results on the stochastic ranking processes has been suc- 
cessfully applied to practical data, such as ranking data of books at an online bookstore 
Amazon.co.jp [HI [13] and list of subject titles at a collected bulletin board 2ch.net jT3] . 

One may wonder why such a simple rule as the move-to-front rule could be observed 
in actual social activities. An explanation is that the ranking numbers on the web (such 
as those representing the books, in the case of online bookstores) usually seek to align 
the web pages in the order of current popularity of the pages. A social impact of the 
development of web-based activities is that it has become possible to catalog a huge 
amount of unpopular items [1]. In fact, a majority of books catalogued on an online 
bookstore are sold less than one copy a month. For such books, any reasonable order 
reflecting the current popularity would be equal to the order of the time of most recent 
sales, because the second recent sale of such book would be long ago, hence would not 
reflect current popularity. Thus the move-to-front rule will provide a simple but universal 
model in the rankings on the web. 

A ranking of a book at Amazon.co.jp jumps close to top of the ranking whenever the 
book is sold at Amazon.co.jp [H], and a subject title in the web page for the list of 2ch.net 
jumps to the top whenever a comment (a 'response') concerning the subject is written 
[T3] . Ordering a book and responding to a subject are social activities which naturally 
are expected to contain day-night difference in the intensity. 

Explicit time dependence, reflecting day-night difference of social activities, are ob- 
served in actual data. Let us regard such time dependence as the non-uniformity of 
intensity measures p-^'*. Pj-^'' are usually unknown quantities to be determined statisti- 
cally from observed data. We then have to consider both particle dependence and time 
dependence in the statistical analysis of the practical data. The assumption of common 
time dependence (11040 developed in Section [5] provides a simple way to take day-night- 
difference of social activity into account, in applying the stochastic ranking process with 
inhomogeneous intensity. 

A.l Factorization of day-night social activity difference. 

In [H], a data taken during the period of about 3 months at Amazon.co.jp is used to 
statistically obtain A, based on (llOSp . The data was taken manually in the year 2007, at 



34 



21:00 each day. We can show that in the case of common time dependence assumption 
fll04p . we can 'factorize' periodic time dependence of a, and that the use of fllOSp in [HI [13] 
is justified in obtaining A from data with periodic time dependence. In fact, assume that 
there exists a positive constant T such that 

(133) fi(t + T) = a(t), t ^ 0. 

We may normahze w^^^^s in f ll04p so that 

1 

(134) - / d{u) du = 1 

T Jo 

holds. Then f ll33p and ( I134p imply j^^'^ {a{u) — 1) du = 0, so that 

(135) Ap{t) := A{t) -t= [ (a(n) - 1) du 

Jo 

is a periodic function with period T, and fllOGp is 

(136) yc{t) = 1 - / e-"'(*+^''(*» X{dw). 

Jr+ 

If we collect data at each fixed time of the day, at tn = tQ + nT, n = 0,1,2, . . ., then 01361) 
implies 

(137) yc{tn) = l- [ e-"'("^+*«+^-(*»»A(c/w). 

Jr+ 

Hence the effect of day-night difference in d is absorbed in the translation of origin of 
time to I— i- to + ^p(to); and the use of formula fllOSp for the constant intensity is justified. 

A consideration of this subsection is of practical use when one has a data much longer 
than 24 hours, as in the case of [14] . 

A. 2 Time change according to intensity measure. 

In [13], a data of list of subject ('thread') titles at a collected bulletin board 2ch.net is 
statistically analyzed using stochastic ranking process. In [13] the data was collected from 
a short period in the daytime, and the problem of day-night activity difference was not 
serious, hence a fit to the formula fllOSp for the constant jump rate (homogeneous intensity) 
was possible [13] . However, to study data of longer periods for sharper statistical results, 
effects of day-night activity difference need to be taken into account. 

In applying f ll04p to the obtained data to extract time dependence (day-night differ- 
ence), we need to estimate the function d in f ll04p or A in f ll07p . This is accomplished by 
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making use of f lTTOj) and (fTT2D . In the case of 2ch.net [13], iV in f lTTOj) or ffTT2|) is about 
700, and since full records of transaction are accessible at 2ch.net, it is possible to put n 
in f lll2p equal to and count all the threads' jumps. In the case of Amazon.co.jp, is 
of order million, and n = N is unrealistic. Even in such cases, if we observe sufficiently 
large number of books (n ^ 1), we can apply the idea introduced here. 

Note that the series Z{N) are approaching their asymptotics in fll30p rather slowly 
for the Pareto distribution. Therefore in practical application of (IllSp with the Pareto 
distribution for A, if one takes A^ = O(IO^) as in 2ch.net [13], one should avoid using the 
asymptotic formula in the right hand side of (11301) . and calculate the finite sums (11161) or 

dniD. 

We announce that we actually collected a 24 hours data of size = 70140 from 
2ch.net, and performed a statistical fit of the data to (11321) . with A^ = 697, and obtained 
b = 0.872 ±0.002. (The error is 90% confidence level. See [1^ for details.) Apparently, we 
have a good single parameter fit to the data, which suggests that the practical assumption 
(11041) is good. Details may be reported elsewhere. 

We note that in [13], a value of 6 = 0.6145 was obtained for 2ch.net (with different set 
of data). This is much smaller than the present result. The data used in [13] was small in 
size, because the data was collected manually in those times, and also, to avoid influence 
of day-night difference in the total activity, the data was for a short time period in [13] . 
so that the result in [13] is less reliable compared to the present result. 

We also note that we have 6 < 1, consistently with previous observation [H] for 
Amazon.co.jp, where we obtained b = 0.809. This shows that, as in Amazon.co.jp, the 
popularity of subjects is concentrated to a relatively small number of threads in 2ch.net. 
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